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SUMMARY 


A number of the most promising integral methods for solving approxi- 
mately the compressible -laminar -boundary -layer equations are investigated 
in order to determine a computationally convenient and sufficiently 
accurate method of calculating boundary-layer characteristics. The chief 
methods considered are: (a) The one -parameter Karman-Pohlhausen method, 

with three different assumptions for the velocity profiles, and (b) the 
two -parameter method, first applied by Sutton, with two different assump- 
tions for the velocity profiles. After the methods are explicitly 
described in general terms for the case of zero pressure gradient and 
for the case of a pressure gradient in the direction of flow with zero 
heat transfer, they are applied to the calculation of the compressible 
laminar boundary layer over a surface with zero pressure gradient, with 
and without heat transfer at the surface, for the purpose of establishing 
the accuracy of the methods. Comparison of the results is made with those 
of known exact solutions for skin-friction and heat-transfer coefficients, 
velocity profiles, velocity derivatives, and especially laminar-boundary- 
layer stability^ From this comparison it is found that the Karman- 
Pohlhausen method with a sixth-degree polynomial as the velocity profile 
is the most suitable for many practical purposes. 


INTRODUCTION 


It is well-known that the differential equations of two-dimensional 
compressible-laminar -boundary-layer flow are difficult to solve exactly. 
Stewartson (reference l) and Illingworth (reference 2 ) have recently shown 
that if the Prandtl number is unity and the viscosity coefficient is 
proportional to the temperature, then the equations for the compressible 
heat-insulated boundary layer with a given pressure gradient can be trans- 
formed into the equations for an incompressible boundary layer with a 
different pressure gradient j however, this principle appears at present 
tedious to apply in practice. 
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The most frequently used and most fruitful methods of solving the 
boundary -layer equations approximately are the integral methods, in 
■which the partial differential equations are integrated over the boundary- 
layer thickness, and are hence satisfied only "in the average." By 
assuming definite forms for the velocity profiles as functions of the 
normal distance, ordinary differential equations are obtained, with 
distance along the surface as the independent variable. Any integral 
method may be regarded as either of two types: (a) The single -integral 

type, in which the partial differential equations are integrated once 
across the boundary-layer thickness, and the profiles contain a single 
parameter to be deter min ed by the resulting ordinary differential equa- 
tion; (b) the multiple -integral type, in which several (say m) integral 
equations are used, and the assumed velocity profiles contain m param- 
eters to be determined by the m resulting ordinary differential equations. 

The best-known integral method is the Karman-Pohlhausen, which is 
of type (a) with fourth-degree profiles. This method has been found 
quite useful for incompressible flow (cf., e.g., Dryden, reference 3 ) 
but it has two important limitations. It fails to predict the separation 
point accurately in an adverse pressure gradient and it often does not 
give sufficiently accurate results for stability calculations based on 
criteria developed by Lin and Lees (references k and 5)*^ For these 
purposes refinements in the usual Karman-Pohlhausen method must be made. 

In the methods of type (a) the refinements usually consist of 
assuming types of profiles which satisfy more boundary conditions than 
the fourth-degree profiles. Schlichting and Ulrich (reference J)j for 
example, have used sixth-degree profiles for incompressible flow, satis- 
fying an additional boundary condition at the wail and also at the outer 
edge of the boundary layer. Satisfactory results were obtained, except 
for flow in the vicinity of a .stagnation point. Weil (reference 8) has 
recently applied this method to compressible flow with zero heat transfer. 
However, no investigation was made here of the expected accuracy of the 
results, although it was pointed out by Weil that the use of sixth-degree 
profiles is expected to yield satisfactory results for stability calcu- 
lations since these calculations involve first and second derivatives of 
the velocity, and the velocity profiles are made to satisfy additional 
conditions involving the rate of change of second derivatives (viz, third 
derivatives) . 

Timman (reference 9 ) bas suggested the use of exponential profiles 
based on exact profiles derived from the solutions of Von Karman and 


-*-Lees (reference 6) has recently applied the Karman-Pohlhausen 
method to an investigation of the stability of compressible laminar 
boundary layers with favorable pressure gradients, but this investigation 
was considered by its author to be essentially qualitative. 



NACA tn 2655 


3 


Millikan (reference 10). Although the calculations are thereby made 
more tedious, it is believed that no significant improvement over poly- 
nomial profiles should be expected, since the latter type of profile 
usually satisfies a fairly large number of conditions at the outer 
boundary-layer edge anyway, and since primary interest usually lies in 
the region near the wall. Moreover, Yuan (reference 11) and Lew (refer r 
ence 12) have also used exponential profiles with no evident improve- 
ment in accuracy. In cases of adverse pressure gradients. Tinman has 
suggested the use of a special profile satisfying an additional condi- 
tion involving the fourth derivative of velocity at the separation point. 
Considerable improvement in the accuracy of prediction of the separation 
point was thereby obtained (reference 9)* 

Loitsianskii (reference 13) has suggested a modification of the 
Karman-Pohlhausen method based on multiplying the momentum equation by 
a small variation of velocity and then integrating across the boundary- 
layer thickness. A velocity profile with a single undetermined param- 
eter is, as usual, assumed. The skin friction is subsequently calculated 
by means of the Karman momentum equation. The method was applied in 
reference 13 for several cases of incompressible flow, but the results 
did not seem to indicate superiority of this method over the usual Karman- 
Pohlhausen method with a fourth-degree profile. 

In the integral methods of type (b) the Karman-Pohlhausen method 
is extended by deriving more than one integral equation. This can be 
done by multiplying the momentum partial differential equation by a 
series of different factors, and then by integrating the resulting equa- 
tions over the boundary-layer thickness. The factors which have usually 
been chosen are integral powers of either the velocity (Leibenson, refer- 
ence 14/ Golubev, reference 15} Sutton, reference 1 6 } and Wieghardt, 
reference 17) or the normal distance (Whitehead, reference 18). A further 
possibility, suggested and applied by Whitehead, is successive integra- 
tion of the momentum equation. The integral equations obtained by these 
procedures are to some extent analogous to those which would be obtained 
by the method of moments, and an infinite number of such equations would 
be equivalent to the original partial differential equation. Because of 
the elaborate nature of the calculations required in such a procedure, 
however, only the first two of the infinite set of equations have usually 
been considered. 

The multiple -integral methods have thus far been developed and 
applied only for incompressible flow. In the present investigation only 
the use' of powers of velocity as factors will be considered. In the 
application made by Sutton (reference 16) a fourth-degree velocity pro- 
file was assumed with two undetermined parameters. However, one of the 
boundary conditions at the wall ordinarily satisfied in the Karman- 
Pohlhausen method was not satisfied. Wieghardt (reference 17) has also 
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used this two -parameter method, hut he has assumed eleventh-degree 
velocity profiles satisfying additional conditions at the wall (in fact, 
the same as those satisfied by the sixth-degree profiles of references 7 
and 8) and at the outer boundary -layer edge. Although results of com- 
paratively high accuracy can thereby be obtained, Wieghardt * s method 
can be quite tedious in practice. This, in fact, is one of the general 
disadvantages of the multiple -integral methods. 

The general case of heat transfer in a compressible laminar boundary 
layer with a pressure gradient is complicated by the fact that in this 
case there is no known solution of the energy partial differential equa- 
tion giving the temperature explicitly as a function of the velocity. 

One means of treating such a case is by transforming the energy, as well 
as the momentum, equation into a differential-integral equation, and 
assuming a profile not only for the velocity but also for the temperature. 
Kalikhman (reference 19) has investigated this case by this means using, 
analogously to the ordinary Karman-Pohlhausen method, fourth-degree pro- 
files for both velocity and stagnation enthalpy. Although important 
useful results are thus obtained, their accuracy is subject to the limi- 
tations of the Karman-Pohlhausen method previously discussed. 

The aim of the present study is to investigate the practical feasi- 
bility of the most promising integral methods (single-integral and double- 
integral methods with various polynomial velocity profiles) from the point 
of view of simultaneous accuracy and ease of computation. The approach 
here is primarily a posteriori. The implications of each of the methods 
considered are developed for only the simplest case, that of a surface 
with zero axial pressure gradient in a subsonic and in a supersonic 
stream, and these are compared with the corresponding implications of an 
exact solution. In particular, results for skin friction and heat trans- 
fer at the surface , velocity profiles, velocity derivatives, and laminar 
stability based on the work of Lees (reference 5) have been considered. 

The comparison of stability is particularly critical, since it is here 
that the largest errors are incurred in the approximate methods, and 
that the greatest differences among the results of the various methods 
appear. 

The following points are among the distinguishing features of this 
study: (a) The two-parameter method involving two differential-integral 

equations is developed and applied to compressible flow, and (b) the 
implications of the various methods are directly compared with respect 
to stability criteriaj this, moreover, is the mhin basis here for judging 
the relative accuracy of the methods for practical purposes. 

Although it does not necessarily follow, without further investi- 
gation, that the conclusions drawn from the analysis of the solutions 
for the flow with zero axial pressure gradient will be valid, without 
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modification, for flows with axi.al pressure gradients, it is "believed 
that the comparisons and resulting conclusions’ should nevertheless serve 
as an indication of the relative merits of the various methods in more 
general cases. 

It may he noted that Mangier* s tranformation (reference 20), which 
brings the boundary-layer equations of axially symmetric flow into the 
form of the two- dime nsional equations, further extends the usefulness 
of two -dimensional -flow solutions. Moreover, axially symmetric flows 
of constant pressure thereby lead to equations analogous to those of 
two -dimensional flow with zero pressure gradient. 

This work, carried out at the Polytechnic Institute of Brooklyn 
Aeronautical Laboratories, was sponsored by and conducted with the 
financial assistance of the National Advisory Committee for Aeronautics. 
The authors wish to acknowledge the helpful discussions of this research 
with Professor R. Paul Harrington and the valuable contribution of 
Dr. H. G. Lew to the formulation of the program in its early stages. 
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SYMBOLS 

coefficients appearing in velocity profile 

11 T 

factor of proportionality in equation lb = C— 

M-oo ^00 

coefficient of average skin-friction drag for surface of 
length x 

specific heat at constant pressure 
specific heat at constant volume 
integrals defined in equations (22) 

local heat-transfer coefficient 
thermal conductivity of fluid 
Mach number 

local Nusselt number (hx/k^) 
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Prandtl number of* gas ^ pCp /k^ 


pressure 

gas constant ^Cp - c^ 

critical Reynolds number 

Reynolds number based on length x 

Reynolds number based on momentum thickness 0 (u^S jv^j 
Sutherland's constant 

transformation variable defined in equation ( 15 ) 
absolute temperature 

velocity at outer edge of boundary layer 
free-stream velocity . 

velocity components in x- and y-directions, respectively 

coordinates parallel and perpendicular to surface, 
respectively 

variable defined in equation (35) 
mass density 

absolute coefficient of viscosity 
ratio of specific heats 

kinematic viscosity (p/p) 

thickness of boundary -layer velocity profile in xy-plane 
thickness of boundary-layer velocity profile in xt-plane 

transformation variable (t/5-j.) 
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Subscripts: 

00 free -stream conditions for subsonic flow} uni form conditions 

behind shock for supersonic flow 

1 local conditions at outer edge of boundary layer 

o conditions at surface 

c values at critical coordinate for stability 

BASIC EQUATIONS 


The following equations describe the steady, two-dimensional, laminar - 
boundary-layer flow of a compressible gas along a surface whose radius 
of curvature is large compared with the boundary-layer thickness: 

Momentum equation in x-direction: 

pu(du/3x) + pv(du/dy) = -(dp/Sx) + (b/dy) [p(^u/Sy)j (l) 

Momentum equation in y-direction: 

$p/dy = 0 (2) 


Equation of continuity: 


(V&0(pu) + (d/dy)( P v ) = 0 (3) 


Equation of state: 


p = pRT 


(4) 


Energy equation: 

puc p (dT/ctx) + pvc p (dT/dy) = u(dp/dx) + (d/<fy)[k(dT/dy)] + p(du/dy) 2 (5) 
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As a consequence of equation (2) the static pressure p at any 
axial position x in the houndary layer is equal to the corresponding 
static pressure in the potential flow, in ■which conditions are assumed 
to he known. 

From the momentum (Bernoulli) equation in potential flow, the 
pressure p can he expressed in terms of the potential -flow velocity 
distribution, thus: 


dp/cix = dp/dx 


- -P1U1V 


(6) 


where the prime denotes differentiation with respect to x, and the sub- 
script 1, conditions at the outer edge of the boundary layer. As is 
usual in aeronautical problems, the coefficients of viscosity and heat 
conductivity are taken to be known functions of the temperature only, 
while the coefficient of specific heat at constant pressure will be 
assumed constant. 

The treatment of these equations is greatly simplified if use is 
made of well-known algebraic velocity-temperature relations which are 
exactly valid for a Prandtl number of unity. These relations are for 
the case of zero heat transfer with or without an axial pressure gradient 



(7) 


and for the case of zero axial pressure gradient with heat transfer 



+ CpT = A + Bu 


( 8 ) 


where A and B are arbitrary constants determined from wall and free- 
stream conditions. With the use of these conditions, equation (7) becomes 


1 



2 


t/t x = 1 + (7 - 1)(mi 2 /2) 


(9) 
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while equation. (8) becomes 


T/Ti = (To/SJl) 



1 + (7 - 



(7 - 1 ) (ml 2 / 2 ) («/°i) 2 



( 10 ) 


The use of equations ( 9 ) and (10) which axe valid only for a Prandtl 
number of 1 in lieu of the more general energy equation ( 5 ) is Justified 
primarily by the simplifications introduced in the analysis. However, 
it has been found that momentum boundary -layer characteristics such as 
skin friction are not significantly affected by a change in Prandtl 
number from unity to the actual value for air, which is between O 065 
and O. 76 . Heat-transfer coefficients are affected approximately 10 per- 
cent by such a change. 

Before discussing the approximate solution of the basic equations 
by means of integral methods, it is convenient to present several sub- 
sidiary relations for later use. Since the pressure within the boundary 
layer at any given value of x is a constant, equation (^) gives 


P/Pl 0 T l/ T 


( 11 ) 


The most accurate representation of the variation of the coefficient of 
viscosity of gases with temperature is usually considered to be Sutherland’s 
formula, namely. 





( 12 ) 


where S is a constant which for air is 216° R and where the sub- 
script 00 denotes reference conditions taken in the undisturbed free 
stream for subsonic flow and in the uniform, flow behind any shock wave 
associated with the leading edge for supersonic flow. Although it is 
possible in principle to use equation (12) in the methods employed here, 
it is inconvenient to do so since certain integrals must then be evaluated 
numerically. It is thus desirable that some acceptable simplifying 
approximation to equation (12) be made. In reference (21) Chapman and 
Rubesin suggested that the mathematical advantages of assuming p and 
T to be linearly related can be obtained along with satisfactory accuracy 
if it is assumed that 
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i*M. - C ( T / T .) 


( 13 ) 


•where C is a factor which is chosen so that the Sutherland viscosity - 
te literature relation is exactly satisfied at the wall temperature T 0 . 
Thus, 

c - (V T -) 1/a ( T « + s )/( t o + s ) 

Equation ( 13 ) will he used throughout this report. 

It is convenient in the method used here to apply the Dorodnitzyn 
transformation. Thus a new variable t is defined so that for a given 
value of x the physical variable y is given by the equation 


y 



dt 


(15) 


Correspondingly the thickness 8^ of the boundary layer in the xt -plane 
is defined as the value of t when y = 5 , and therefore 


8 



( 16 ) 


With the additional relations given by equations (ll), (13) * 

(l4), and ( 15 ) , equation (l) can be recast into an infinite set of 
integral-differential equations. This can be done by multiplying equa- 
tion (l) by u 11 dy(n = 1,2,3* . . .) and by integrating from y = 0 to 
y = 8 with the boundary conditions that at y = 0 


u = v = 0 


U = U n 


Su/dy = 0 




(17) 


and at y = 8 
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After transforming to the xt -plane, nondimensionalizing, and introducing 
r = t/S^, the following set of equations is obtained: 


(n + lJpjU,; 


n+2 


-1 


5 t (d/dx) 


pi u i” 2B t - W a 1) 


m+li jl 


K/ te )( 6 t 2 K) f^[ u l u i - ( T / T a) ( u / u i) n 


dr = 


^ Jl( p l U l T ») _1 J 0 ( u / U i) n ' 1 L(^ T) (u/u i )_ 


C^P^T^" 1 (8/dr) (u/uj 


dT n > 1 


n = 0 


( 18 ) 


The application of the infinite set of equations represented by 
equations (l 8 ) will now be discussed. By means of equations ( 9 ) and (10) 
the temperature can be expressed in terms of the velocity component u. 
Wow if it were assumed that 


u/Uq = a^x^J (19) 

J=0 


then several of the a^ coefficients could be determined in terms of 

the other coefficients so that certain boundary conditions, at least 
those in the xt-plane corresponding to those given by equations ( 17 ) in 
the xy-plane, would be satisfied. The rest of the coefficients and 5-j. 

could be calculated from the infinite set of first-order ordinary differ- 
ential equations represented by equations (l 8 ). Such a solution would 
in principle be exact and would be equivalent to a solution of the 
original partial differential equations. 

In actuality, however, since the attainment of such a solution 
would be quite difficult, it will usually be considered sufficient for 
practical purposes to obtain approximate solutions to the infinite, set 
by taking only the first few differential equations of this set and 
correspondingly only the first few terms of the power series in equa- 
tion (19). Indeed, much information of engineering importance for the 
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flat-plate case has been obtained in the past by using only the equation 
corresponding to n = 0, that is, to the Von Karman momentum integral 
(e.g., references 3> U, and 12). Additional accuracy may be obtainable 
by the use of a second equation, that is, by the use of equations (l8) 
for n = 0 and n = 1. Therefore, these two equations will now be 
explicitly written. 

For n = 0: 

( F l/ 2 )( s t 2 ) ’ + |V + Pl u l 2 )' - = 

( F 3P- CT l)/(Px t, l T -) f 20 ) 

and for n = 1: 

K/ 4 )( s t 2 )' + [(V/ 2 ) + ( F u/ 2 ) ( lo Se Pl«l 3 )' - F 5 (loge u x)] H 2 " 
(( F 6P» c )/(Pl°l)]( T l/ T .) I 21 ) 

where the prime denotes differentiation with respect to x and where 

p i - § - ( u /°i)] ar 

* 2 =/ 0 [(wO - ( t / t 4] 4t 

*3 - [(5/3t) (uK)] 0 

= J_ ( u / u l) [x - ( u / u x ) 2 4 t 
f 5 = f (W u x)R - ( T / T x)] 
f 6 = [ |(p/3t) (u/Xxjj 2 dr 
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In these Fj integrals either equation (9) or (10), as the case may 
he, is used to express T/Tq. 

In the actual application of either or both equations (20) and (21), 
u/Uq is assumed as some function of t. For convenience a polynomial 

in t is here assumed, that is 


u /° 1 




(23) 


where N is the degree of the assumed polynomial containing (N + 1) of 
the a j 1 s . 

The a j * s are determined from the boundary conditions and from the 

differential equations (20) and (21). The boundary conditions in the 
xt-plane corresponding to those given by equations (17) must be satisfied. 
Furthermore, additional accuracy can be obtained if the a^'s are selected 

so that the approximate solution given by this method has the same value 
of lower derivatives at y = 0 and y = 5 as an exact solution to the 
partial differential equations would have. These values may be obtained 
by differentiating equation (l) one or more times with respect to y 
or t. For completeness the useful boundary conditions in the x,t coordi- 
nates will now be listed for flows having zero axial pressure gradient 
with or without heat transfer or having an axial pressure gradient but 
no heat transfer: 


At t = 0: u/Uq = 0 

^~ C ( aS / S t 2 ) ( u / U l) = -Pl u x’ 6 t 2 
(bVcH-3) (u ju-^j = 0 

At t = 5 t : u/Uq = 1 

(u/Uqj = 0(m = 1,2,3, • . .) 


( 24 ) 


It may be pointed out that although additional conditions at y = 0 
can be derived they become impractical ly cumbersome to apply. 
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The method involv ing the simultaneous solution of equations (20.) 
and (21) will be referred to as the "two -parameter method," while that 
of using only equation (20) wl.ll be the "one -parameter method." In the 
two -parameter method one coefficient a^ in addition to 6^. is deter- 
mined from the two differential equations, the rest from the boundary 
conditions. In the one-parameter method S-j. is determined from the 

solution of equation (20) and all of the a^ coefficients are deter- 
mined from boundary conditions. 

In carrying out the integration of these differential equations 
graphical or numerical methods are in general required. The arbitrary 
constants are usually determined so that the -unknowns are either well- 
behaved or have definite prescribed values at some value of x. 

Before applying the methods indicated here to the case of an axial 
pressure gradient, it is considered de sizable to investigate critically 
the accuracy obtained by using several types of profiles in equation (20) 
and in both equations (20) and (21) for the compressible flow with no 
axial pressure gradient. An exact solution for this case has been given 
by Chapman and Rubesin (reference 21). The accuracy of the results will 
be measured by two criteria, namely, by the accuracy of the skin-friction 
and heat-transfer values and by the critical Reynolds number for the 
stability of the laminar boundary layer. 

The application of these integral methods to the flat plate and 
the results for skin friction and heat transfer will now be discussed 
in detail. 


BOUNDARY LAYER WITH ZERO AXIAL PRESSURE GRADIENT 


The one -parameter and two -parameter methods described in the pre- 
ceding section are here explicitly applied to the compressible flow over 
a surface with zero axial pressure gradient. Various profiles will be 
assumed and the results compared with those obtainable by a mathematically 
exact solution for this case. 


Solution of Equations 

For the case of zero axial pres sore gradient U^ = U and equa- 
tions (20) and (21) become 


(?t 2 )’ * (2V^ 3 )/(u,J?i) 


( 25 ) 
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and 

( 6 t 2 )' - (K cf 6)/('”\) ( 26) 


where the F^ integrals are in this case constants, and not functions 
of x. 

In the two-parameter method equations (25) and (2 6) are solved for 
and the additional a^ coefficient which has not teen determined 

from boundary conditions. With the initial condition S t = 0 at x = 0 
these unknowns are determined by the equations 

2Fg/F 4 = F3/F-L (27) 

and 

s t/* = N /( 2F 3/ F l) fK) (28) 

where = U^x/v^. Equation (27) leads to a quartic equation in the 
unknown coefficient a^, three roots of which must be rejected on physical 
grounds. Once the constant value of a^ is known, F3 and F^ can 
be evaluated and substituted into equation (28) to complete the solution. 

In the one-parameter method the F^ integrals are completely deter- 

d 

mined once the velocity profile is chosen so as to satisfy the selected 
boundary conditions. Thus, only equation (28) is used to determine the 
single unknown 5^ . 


Skin Friction and Heat Transfer 


From the solutions outlined in the preceding section, the skin- 
friction and heat-transfer characteristics are readily determined. The 
coefficient of average skin-friction drag C f is defined by the equation 


C 


f 




(29) 
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By use of equation (28) this can he written as 

C f = 2 ^^3 (c/i^) 1 / 2 (30) 

Similarly a local heat-transfer coefficient h can he defined as 

h = -(k dT/dy) 0 |(T o - T e ) (31) 

where T 0 is the equilibrium wall temperature for no heat transfer and, 
for Npj. = 1, is given hy the equation 



By the use of equation (28) the equation for h is found to he 

%u = 

= (^3^/2) 1/2 (33) 

Hence, hy comparison of equations (30) and ( 33) ^ 

%u = c Al h (34) 

Velocity Profiles and Results for Skin Friction and Heat Transfer 

Several velocity profiles were assumed and the friction and heat 
transfer given hy these profiles were calculated. In the two-parameter 
method fourth- and fifth-degree polynomials were used, while in the one- 
parameter method fourth-, fifth-, and sixth-degree polynomials were 
chosen. 

The results obtained here are compared mainly with those of the 
exact analysis of Chapman and Rubesin (reference 21) which treats a 
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general Prandtl number and variable surface temperature. The exact 
velocity profiles of reference 21 are derived from the well-known Blasius 
differential equation, solutions of which are tabulated, for example, 
in reference 22. With the assumptions Np r = 1 and constant surface 

temperature for which equations (9) and (10 ) are valid, exact temperature 
profiles are thus obtainable from the exact velocity profiles of refer- 
ence 21. (These profiles, when expressed in terms of the Blasius vari- 
able Tj (cf. equation lb-, reference 21), are independent of the surface 
temperature and the Prandtl number. ) 

The variation of skin friction with Mach number for an insulated 
surface is represented in figure 1. Included for comparison are results 
already calculated in reference 21 for Np r = 0.72, as well as exact 

results for Np r = 1. 

Table I gives values of the skin-friction coefficient and Nusselt 
number obtained by the one- and two-parameter methods with the various 
velocity profiles. Values derived from the exact analysis of refer- 
ence 21 for Np r = 1 and 0.72 axe listed for comparison. The free- 

stream temperature and the wall temperature, which for an insulated 
surface is a function of the free -stream Mach number, are contained 
explicitly in the term C. Typical values of Vc for an insulated 
surface are presented in table II in connection with the calculation of 
the skin-friction drag coefficient. 

Table I shows that the values obtained for skin friction and heat 
transfer by the two-parameter method with either of the two profiles 
used and by the one-parameter method with a sixth-degree profile differ 
from the exact values by less than 0.6 percent. The values obtained by 
the one-parameter method with the fourth-degree profiles are in error 
by approximately 3 percent. Thus, all methods used here give sufficiently 
accurate results for engineering purposes for skin-friction and heat- 
transfer coefficients. 

Since in both the exact solutions (reference 21) and the solutions 
obtained here by the integral methods, the effect of Mach number in the 
case of an insulated plate with Np r =1 is given by the same factor, 

namely \fc, it follows that (for a linear temperature-viscosity relation) 
these integral methods lead to the exact effect of Mach number on Cp 

and I% u . It may be noted that the Mach number effect depends solely on 

the assumed temperature -vi sco sity relation. The effect of the latter 
can be seen from figure 1, where the results of the various integral 
methods and of reference 21, all based on a linear temperature -viscosity 
relation, are shown along with Crocco’s results which are based on the 
Sutherland formula. 
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From figure 1 it can "be seen, incidentally, that the value of the 
Prandtl number has a smal l effect on the variation of Cf with Mach 

number in the case of zero heat transfer. For example, in the range 
M = 0 to M = 5, with T m = 6k8° R, C f] /B^ varies froni 1.328 to O .983 
for Rp r = 0.72, and from 1.328 to 0.95^ for Hp r = 1. It can be seen 

from table I, moreover, that the Prandtl number has a somewhat larger 
effect on the value of Nusselt number at all Mach numbers, the change 
from Np r =1 to Np r = 0.72 being about 10 percent. 


Comparison of Velocity Profiles and Derivatives in Physical Plane 

For estimating the stability of the laminar- boundary layer, the 
accuracy of the first and second derivatives of the velocity profile, 
as well as that of the profile itself, are of importance. In this sec- 
tion the profiles and derivatives obtained by the integral methods under 
discussion are con^ared with the exact Blasius solution of references 21 
and 22. The comparison can be conveniently made in the physical xy-plane 
by introduction of the dimensionless variable 


6 = (y/2x)(R x /c) 1 /2 


(35) 


The derivatives with respect to e are shown in figures 3 and- 

It is seen from figures 2, 3 } and 4 that, although the velocity 
profiles obtained by the integral methods agree well with the exact 
profile, the derivatives, as might be expected, do not correspond so 
closely. Previous comparison of skin friction and heat transfer, which 
depend upon the first derivatives at the surface, indicates that each 
of the integral methods gives good results for these characteristics. 
Hence each of the integral methods may be considered to predict the 
value of the first derivative at the surface with sufficient accuracy. 
Moreover, except for the fourth-degree profiles, the first and second 
derivatives given by a ll the approximate solutions are in good quali- 
tative agreement with the exact results throughout the boundary-layer 
thickness. Examination of figures 3 and ^ for first and second deriva- 
tives over the boundary -layer thickness indicates, however, that the 
one-parameter method with a sixth-degree velocity profile gives on the 
whole the most satisfactory results. It is to be anticipated that this 
method will consequently also give satisfactory results for stability 
calculations . 
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DETERMINATION OF CRITICAL REYNOLDS NUMBER FOR STABILITY 


The experimental work of Schubauer and Skramstad (reference 23) has 
clearly established that the transition from a laminar to a turbulent 
boundary layer is due to an instability of the laminar layer, if the 
turbulence in the free stream is low, and if the surface of the body 
is aerodynamically smooth and has a large radius of curvature. Thus 
the laminar-boundary-layer stability theory, which has been developed 
for incompressible flow by several investigators over the past 25 years 
and which has been recently extended by Lin and Lees (references 4 and 5) 
to the practically interesting case of compressible flow, may be used 
to predict the local Reynolds number at which instability of the laminar 
layer will first occur, or above which self -propagated disturbances are 
not damped out. Transition to turbulent flow will take place downstream 
of the point corresponding to this local Reynolds number, which is termed 
the critical Reynolds number j the exact distance downstream cannot be 
predicted by the present stability theory, which is based on small- 
perturbation methods, but seems to be dependent on the value of the 
critical Reynolds number, on the magnitude of the small but possibly 
finite turbulence in the free stream, on the surface roughness, and on 
the potential -flow pressure gradient. 

The exact calculation of the critical Reynolds number for a given 
velocity profile is tedious. However, the approximate stability rules 
of Lin and Lees (references 4 and 5) permit a rapid determination of 
this critical value with a minimum of labor. These rules have been 
shown by the work of Lin and Lees and that of Hahneman, Freeman, and 
Finston (reference 2k) to give reliable results and have, therefore, 
been used extensively. It might be further mentioned that because of 
approximations inherent in the laminar-boundary-layer stability theory, 
and because of the undetermined relation between boundary-layer neutral 
stability and transition, there seems to be no practical justification 
for obtaining greater accuracy than that given by the approximate sta- 
bility rules. 

The results of stability calculations are sensitive to the accuracy 
of the profiles. Thus the critical Reynolds numbers based on approximate 
solutions to the boundary-layer equations have not in general been in 
good agreement with those based on exact solutions. The reason for this 
discrepancy is clear when it is considered that the stability calculation 
depends strongly on first and second derivatives of the velocity profile. 
Approximate solutions may be expected to give fairly accurate velocity 
profiles, but unless special care is exercised the first and second 
derivatives of the profile throughout the entire boundary-layer thick- 
ness might be quite in error. This has been shown in the previous sec- 
tion where it has been pointed out that the frequently used Karman- 
Fohl hausen method with a fourth-degree velocity profile gives reasonably 
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accurate values for the skin friction, the heat transfer, and the velocity 
profile hut yields quite inaccurate values for the derivatives. It may 
he -expected, therefore, that a crucial test of the accuracy of approxi- 
mate solutions to the laminar -houndary-layer equations would he a com- 
parison of the critical Reynolds numbers predicted hy these solutions 
with those predicted hy exact solutions. 

In this section the critical Reynolds numbers of the profiles 
obtained in the previous section for the zero-heat-transfer case are 
obtained and compared with those of the exact solution of Chapman and 
Rube sin. 

With the prime denoting partial differentiation with respect to 6, 
Lees’ criterion modified for the viscosity-temperature relation used here 
'is given as follows: 


E cr - 2 (^,cr c ) l/2 

- 2 5(”/Wl) o , °(t/ 5 !0 2 (“/%)' 4 {i - m i 2 £ - ("/ u i)^ 2 }' l/2 ( 36 ) 


where the subscript c denotes values at the point 6 = e c at which 
the following relations are satisfied: 


J - -*( T l/ T o) (»/%)o'("/»l) [( u /“l)] " 3 |(f/ T l)( u / D l)'' - (»/“l) '(?/%)] 

(37) 

and 


O. 580 /J = 1 


2 


( U / B l)c' e ( u / u i) 


1 


( 38 ) 


It will be noted from equation (37) that the first and second deriva- 
tives of the velocity and temperature profiles in the critical region 
around e c are influential in the determination of R cr . The procedure 

for obtaining R cr is to calculate J from equation (37) for various 
values of e, and to find the values of e = e c for which equation ( 38 ) 
is satisfied. After the values of (u/u -j_) c and (T/T-jJ c are determined 
at this point, R cr and Rj. cr follow from equation ( 36 ). 
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For comparison with previous results it is convenient to calculate 
the critical Reynolds number Rg cr "based on momentum thickness. Since 

Q is defined by 


0 = f (p/pi)( u / u i)[ 1 - ( u / u i)] d y 


(39) 


it is found, for the integral methods, that 

Be.cr - ( F l F 3/ 2 ) l/% cr 
and for the profiles of reference (21) 

R 0,cr = 0*332R cr 


m 


(hi) 


Calculations of the minimum critical Reynolds numbers have been 
made for the flow over a flat insulated surface at free-stream Mach 
numbers of 0, 1, and 2. The results are tabulated in table III and 
plotted in figure 5. Primarily, comparison of the integral methods 
presented here should be made with the results derived from the exact 
profiles of reference 21 since identical temperature-viscosity relations 
are used in both. The results of Lees, given in figure 5 of reference 
are also presented for comparison. Note should be taken that Lees’ 
approximate values are based on still a further approximation to equa- 
tions (36), (37 ) > and (38) presented here. 

Figure 5 indicates that the one -parameter sixth-degree polynomial 
method gives the best general agreement with results of the exact solu- 
tions over the range of Mach numbers studied. The two-parameter, fifth- 
degree method is next best, while in order of decreasing accuracy, the 
one-parameter fifth-degree, two-parameter fourth-degree, and one -parameter 
fourth-degree polynomials follow. The two -parameter methods seem to 
give better general agreement than the one -parameter methods for equal 
degree of the assumed polynomial for the velocity profile. (Because of 
the inaccuracies introduced by the approximation in this stability 
criterion, as well as by the required graphical method involved in its 
application, the derived results are not reliable to more than two signif- 
icant figures. ) 
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The stability analysis given here has thus indicated that integral 
methods of solution of the houndary-layer equations may he used to 
predict the critical Reynolds number provided that either additional 
integral equations or higher-degree polynomials and additional boundary 
conditions with a single integral equation are used. Thus the choice 
of method to he used in the more general cases of pressure gradient with 
and without heat transfer is either the two-parameter method or the one- 
parameter method with higher-degree polynomials, particularly of the 
sixth degree. 

While these two methods lead to results which are of nearly the 
same accuracy, use of the one -parameter method has the advantage of 
simplicity of calculation, which becomes significant in cases involving 
pressure gradient. It is thus concluded that the one-parameter method 
with sixth- or higher-degree polynomials is the most satisfactory and 
promising method for extension to general compressible boundary-layer 
analyses. Sixth-degree polynomials will probably prove to be satis- 
factory in many cases, but use of a seventh-degree profile satisfying 
an additional condition at the separation point (cf. reference 9 ) may 
increase the accuracy of determining the separation point in adverse- 
pressure -gradient cases. Moreover, a sixth-degree profile may not be 
satisfactory for flow near a stagnation point (cf. reference 7 ). 


CONCLUSIONS 


The laminar -boundary-layer equations for compressible flow can be 
converted into one, two, or more integral-differential equations. Approxi- 
mate solutions can be obtained by assuming special forms for the velocity 
profiles satisfying various boundary conditions and containing, in 
practice, either one or two parameters to be determined by these equations. 

From a comparison of the integral methods discussed here for the 
compressible flow over a surface with no axial pressure gradient, the 
following conclusions can be drawn: 

1. All of the methods predict the values of the skin-friction and 
heat-transfer coefficients at the wall as well as the over-all velocity 
profiles with satisfactory accuracy for engineering purposes. The two- 
parameter method and the one -parameter method with a sixth-degree pro- 
file are particularly accurate . 

2. The one-parameter method with a sixth-degree velocity profile 
gives on the whole the most accurate results for velocity profiles and 
their derivatives throughout the boundary-layer region. 
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• 3 . The critical Reynolds number for laminar -boundary-layer sta- 
bility is predicted with qualitative accuracy by all methods. Moreover, 
reasonably good quantitative accuracy is obtained with the two -parameter 
method and with the one-parameter method with a sixth-degree profile. 

4. Because of computational simplicity and equality of accuracy, 
the one -parameter method with sixth- or higher-degree polynomials appears 
to be the most satisfactory and promising method of investigating the 
more general case of laminar compressible boundary layer with axial 
pressure gradient and heat transfer. 


Polytechnic Institute of Brooklyn 

New York, N. Y., November 6, 1950 
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TABLE I.- SKIN -FRICTION DRAG COEFFICIENT 
AND NUSSELT NUMBER 




























TABLE III.- MINIMUM CRITICAL REYNOLDS NUMBERS TOR STABILITY 
IN LAMINAR -BOUNDARY -LAYER PLOW OVER AN INSULATED 


Method 

(1) 


Lees 1 exact 
stability 
calculation 

Lees ' approximate ' 
stability 
criterion 



M« - o 


Rnv» 


no 27,000 


0.4156 


0.499 n9 32,000 
0.945 .517 290 96 21,000 1.97 O.72I 

.862 .47 46o 152 53,000 2.05 .75 

.958 .51 310 102 24,000 2.17 .77; 

00 .829 .46 520 178 68,000 2.05 .73 


Chapman-Rubesin 0.64 . 419 540 180 

2 parameter-4 .56 .370 890 290 


2 parameter -5 
1 parameter -4 


1 -rio-r*orm»+rwr 


.56 .370 890 290 

• 59 *381 780 260 

.52 .345 1210 420 

Ao 4-56 k<V> T k^ 


1 parameter-6 .61 .394 69O 230 


■For example, 2 parameter-4 means two-parameter method with a fourth-degree profile. 


e c 

: G)J 

1.97 

0.724 

2.05 

.75 

2.17 

.773 

2.05 

• 73 

2.20 

7flA 

• 1 w— 

2.07 

.75 


M« = 2 


R-y> 


1480 
960 
630 
1480 
480 
21 960 
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(a) T„ = 6bQ° F absolute. 


C~ R , N== l 


I- 6 


E-4.5 


(b) T ro = 72° F absolute. 


Figure 1.- Variation of skin-friction drag coefficient with free-stream 
Mach number for insulated flat surface. I, one-parameter method; 

II, two-parameter method; 6, degree of polynominals; C-R, 

Chapman-Rube sin (reference 21); C-S, Crocco's calculation based on 
Sutherland's temperature -viscosity relation (cf. reference 21). 
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(a) M ro = 0. 

Figure 2.- Comparison of velocity profiles for flow over insulated flat 

surface . 


ONE- PARAMETER 
SIXTH-DEGREE POLYNOMIAL 
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(a) Two-parameter method; M ra = 0. 


Figure 3 .- Comparison of first derivatives of velocity profiles for flov 

over insulated flat surface. 




(b) One -parameter method; = 0. 


Figure 3.- Continued 
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(c) Two-parameter method; M ro = 2. 


Figure 3 .- Continued. 
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(a) Two-parameter method; M„ = 0. ' 

Figure - Comparison of second derivatives of velocity profiles for flow 

over insulated flat surface. 
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(c) Two-parameter method; = 2. 
Figure 4.- Continued. 
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(d) One -parameter method; M m = 2 
Figure k.- Concluded. 
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1*0 




Figure 5*- Variation of minimum critical Reynolds number with Mach number 
for laminar boundary layer flow over insulated flat surface with 

»Er = 


NACA - Langley Field, Va. 




